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Abstract 

The thermal and phase properties of a multifragmentation model which uses clusters as degrees 
of freedom, are explored as a function of isospin. A good qualitative agreement is found with 
the phase diagram of asymmetric nuclear matter as established by different mean-field models. In 
particular, from the convexity properties of the nuclear entropy, we show that uncharged finite 
nuclei display first and second order liquid-gas-like phase transitions. Different quantities are 
examined to connect the thermal properties of the system to cluster observables. In particular 
we show that fractionation is only a loose indication of phase coexistence. A simple analytical 
formula is proposed and tested to evaluate the symmetry (free) energy from the widths of isotopic 
distributions. Assuming that one may restore the isotopic composition of break-up fragments, it 
is found that some selected isotopic observables can allow to quantitatively access the freeze-out 
symmetry energy in multifragmentation experiments. 

PACS numbers: 25.70.Pq Multifragment emission and correlations, 24. 10. Pa Thermal and statistical models 
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I. INTRODUCTION 



Due to the short range repulsive and finite range attractive character of the nucleon- 
nucleon interaction, nuclear matter is known to exhibit a phase transition similar to the 
liquid-gas transition taking place in real fiuids l2|. 

Many experimental and theoretical efforts have been devoted to this subject js]; in par- 
ticular, most of the recent works focus on asymmetric matter and the effect of isospin as an 
additional degree of freedom . Mean-field based models have demonstrated the 

presence of a first order phase transition in both isospin symmetric and asymmetric nuclear 
matter Sj, the decrease of the critical temperature with increasing isospin asymmetry and, 
in the case o 
gas phases {4, 

In astrophysics, both neutron-star structure and supernova dynamics are infiuenced by 



asymmetric matter, different neutron/proton compositions of the liquid and 



thermal properties of neutron-rich nuclear matter in a large interval of temperatures and 
densities jol. llo|. All these different phenomena involve excited matter at baryon densities 
lower than normal nuclear matter density: this corresponds in the phase diagram to a 
region of instability with respect to phase separation. Information on the phase structure 
and properties of hot and diluted nuclear matter is thus of clear astrophysical relevance. 

From an experimental point of view, the only terrestrial phenomenon that may allow to 
access finite temperature low density properties of neutron rich matter, and in particular pin 
down the density and temperature dependence of the symmetry energy, is given by nuclear 
multifragmentation Indeed, this specific decay channel of nuclei, whose excitation energy 
is of the order of a few MeV/nucleon, has been since a long time tentatively associated to 



the coexistence zone of the nuclear matter phase diagram 

However, the connection between multifragmentation observables and the nuclear matter 
energy-density functional and phase structure is far from being trivial. First, finite nuclei are 
charged while nuclear matter is by definition neutral. Therefore, the presence of this non- 
saturating long range force makes difficult to rely the multifragmentation phenomenology 



to the theoretical studies of nuclear matter 



56!]. Second and even more important, atomic 



nuclei are composed of a very small number of constituents, and their phase properties are 
not trivially linked to the phase structure of nuclear matter. At the thermodynamic limit, 
the coexistence zone is a simple linear superposition of pure liquid and gas phases, and 
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can be deduced from a mean field approach through a Gibbs construction. The situation 
is completely different in finite nuclei, where phase coexistence is revealed by convexity 
anomalies of the entropy surface 



13 



14| . The properties of coexistence cannot be deduced 



from the properties of the pure phases (nuclei and nucleons, respectively), mean- field based 



approaches 
mandatory 



Dadly fail, and a description explicitly accounting for complex clusterization is 



13|. 



A typical example of this ambiguity is given by the fractionation phenomenon, originally 
expounded as a consequence of Gibbs phase equilibrium, later on recognized as a generic 
uster formation 1^, 16, 17 1, and systematically observed in the experimental 



feature o: 
analyses 
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22l |: a clear thermodynamic interpretation of fractionation is not 



trivial in a finite system since we cannot unambiguously associate a given fragment size to 
the "liquid" or "gas" phase. 

To contribute to match the gap between nuclear matter thermodynamic studies and nu- 
clear multifragmentation, in this paper we investigate the isospin dependent phase diagram 
of finite excited nuclei. The study is done in the framework of the Microcanonical Multi- 
fragmentation Model (MMM) {2311, that explicitly considers clusters as degrees of freedom. 
A similar analysis was already presented in Ref. [24i], where Coulomb effects on the phase 
diagram were especially addressed. Here, in order to concentrate on isospin effects, we first 
consider that the Coulomb interaction is switched off. The use of such an idealized neu- 
tral system will additionally allow us to make connections with the expected behavior at 



;he thermodynamic 
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271, 
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imit and to explore the link between isospin observables (isoscaling 



291], isospin fluctuations [18 
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321]) and the low density finite 



temperature symmetry (free) energy of the equation of state. The robustness of measure- 
ments of the symmetry energy in both finite neutral systems and real nuclei under the effect 
of mass and charge conservation and Coulomb is finally addressed. 

II. THEORETICAL FRAMEWORK 



To provide a realistic description of multifragment production, a theory dealing with com- 
plex correlations well beyond the mean field is needed. An exact solution of this problem at 
the microscopic level is provided by classical models (Molecular Dynamics, Lattice Gas) ll| , 
which however completely miss all the specific quantal features of the process. Moreover, no 
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connection is possible within these models between fragment observables and the different 
ingredients of the nuclear energy-density functional, which are explored through multifrag- 
mentation reactions. On the other hand, semiclassical or quantum molecular models (such as 



QMD 



331], AMD 3J]) which reproduce the most important macroscopic nuclear properties 



as density distributions and binding energies, and account for nucleon-nucleon interaction 
and Pauli blocking, still do not offer a definite description of multifragmentation and access 
to the equation of state of excited matter because of the ambiguity of the fragments and 
break-up definition. 

An interesting alternative is given by statistical models which use clusters as degrees of 



freedom 



35| . Such models offer the remarkable advantage that all nuclear bound as well as 



continuum states are naturally accounted for via empirical parameterizations of the clusters 
energies and level densities, allowing a direct comparison with experimental data. 

The price to pay for such a realistic inclusion of nuclear effects, is the underlying hy- 
pothesis that nuclear correlations are entirely exhausted by clusterization, which amounts 
to the implementation of the properties of isolated low excited nuclei for the description 
of break-up fragments. This limitation can in principle be avoided by including effective 



361]. In this case, however, the 



in medium corrections in the fragments energy functional 
fragment energy becomes a free parameter of the theory. 

The distinctive features of the fragmentation process, namely the explosion of an isolated 



nucleus into vacuum, recommends the microcanonical framework [37|, 138|, 139| as the most 
natural choice. The unphysical hypothesis of a sharp fixed freeze-out volume constraint may 
be easily overcome by considering a total spatial extension for the fragmenting system fluc- 
tuating event by event [4^ . Technically, this is realized introducing a A Lagrange parameter 
conjugate of the volume V which alters the statistical weight of a configuration Wc by an 
extra factor, exp(— AV^) [41]. The thermodynamical potential associated to this ensemble is 
Se [X] = S-XV = \nW{E, A), where W{E, A) = / W{E, V) exp{-XV)dV. In addition, for 
a system belonging to the liquid-gas universality class, the exploration of the configuration 
space along constant A paths provides a straightforward method to reveal phase coexistence 
by the back-bendings of the corresponding caloric curves and to finally construct the phase 
diagram |l4]. 



The MMM version [231 of the microcanonical multifragmentation models 



37 



38 



39] 



has been used so far to investigate the thermodynamic properties of charged nuclei with 
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excitation energies between 1 and 15 MeV/nucleon 2J], and will be presently employed for 
the study of isospin effects. MMM provides a Monte-Carlo calculation of the global density 
of states Z, E, P, L, V) of a nuclear system modelized as a non-interacting collection 

of nuclear clusters. The space of observables is given by the baryonic number A, the proton 
number Z, total energy E, total momentum P, total angular momentum L and freeze-out 
volume V. The investigation of all clusters states compatible with conservation laws and 
geometrical restrictions is performed using a Metropolis trajectory in configuration space. 

Break-up fragments are considered as having normal nuclear density Pq, and described by 
a ground-state liquid drop binding energy including surface and symmetry^ terms. This de- 
scription is consistent with a semi classical Thomas- Fermi approximation [42] or hot Hartree- 



Fock 



431], where the effect of temperature is a modified occupation of the single particle 



eigenstates of the mean field hamiltonian. The finite temperature fragment energy func- 
tional in this approach is thus modified respect to the ground state only for the internal 
excitation energy (e) coming from the occupation of continuum states, which are treated 
with a Fermi gas level density parameterization. To avoid double counting of the free par- 
ticles states [44], a high energy cut-off (r=9 MeV) is applied to the level density. 

To allow a comparison with the well known nuclear matter thermodynamics j^, ^ Q] and 
best isolate isospin effects on the fragmentation process, we ignore the long range Coulomb 
interaction and, in order to avoid interference with finite size effects, consider equal size 
systems which differ by the neutron/proton ratio. 



III. PHASE DIAGRAM 



The isospin dependence of the phase diagram for the MMM model is easily spotted in 
the microcanonical "isobar" ensemble, where energy is fixed and volume fluctuations are 
allowed and controlled through a conjugated Legendre intensive. 

Indeed, contrary to ordinary (macroscopic) thermodynamics, the thermal and phase char- 
acteristics of a finite system depend on the statistical ensemble considered. The liquid-gas 
phase transition has a non-zero latent heat and has density as an order parameter, meaning 
that the two associated phases can be distinguished by their different particle and energy 
densities. If a finite system (with a given fixed particle number) exhibits this transition, 
its event distribution at the transition point will show the two peaks corresponding to the 
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two phases if and only if both energy and volume are free to fluctuate, i.e. in the canonical 



isobar ensemble 



1^ 



The general relationship between a distribution in the ensemble characterized by an 
intensive variable 7 associated to the conjugated extensive variable m, and the Boltzmann 
entropy W{m) = exp S{m), 

P^(m) = Z^^ exp (S{m) — 7m) , (1) 

insures then that at the liquid-gas transition point the compressibility is negative in the 
canonical isochore {f3, V) ensemble, and the heat capacity is negative in the microcanonical 
isobar [E, P) ensemble. 

In the multifragmentation transition described by MMM, the low multiplicity ordered 
phase (compound nucleus) and high multiplicity disordered phase (multifragmentation) can 
be distinguished by their energy and volume, just like in regular liquid-gas, and all the above 



considerations apply j24l |. 



Fig. [T] shows some constant A microcanonical caloric curves of a 200-nucleon systems with 
different neutron-proton ratios. The Lagrange parameter A can be associated to a pressure 
through P = AT, where T = [OSe [A] /dE^ ^ is the constant A microcanonical temperature. 
The expected isospin invariance in the absence of Coulomb is confirmed by the fact that 
mirror nuclei ((200,70) vs. (200,130), and (200,50) vs. (200,150)) show an identical ther- 
modynamical behavior, translated into fully superimposable caloric curves. The symmetric 
system shows a broad back-bending, signaling a liquid-gas-like phase transition. Increasing 
the isospin asymmetry, the temperature shows a monotonic decrease, and the back-bending 
width shrinks. The A value A = 3 • 10~^ fm~^, which still corresponds to the coexistence 
region for the Z/A = 0.35 system, appears clearly super-critical for Z/A = 0.25. 

The presence of Coulomb effects in physical nuclear systems breaks the n-p invariance 
and can considerably mask the isospin effects shown by Fig. [1] [24]. In order to disentangle 
these effects, most of the experimental analyses concentrate on different isotopes of the same 
element (eg. ^^^Sn+^^^Sn and ^^'^Sn+^^^Sn collisions). In this case however, the sources have 
different total sizes and a naturally rising question is to what extend finite size effects inter- 
fering with asymmetry effects may blur the signals of the last ones. A quantitative answer is 
offered by Fig. [H where caloric curves corresponding to symmetric systems 50% larger and 
50% smaller than the previously discussed ones are considered. The relative displacement 
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E (MeV/nucleon) 

FIG. 1: (Color online) Microcanonical caloric curves at constant A = (/3P) = 3 • lO^^fm"^ for 
200-nucleon systems and various neutron-proton ratios as indicated in the legend. E is the total 
energy of the system. The magnitude of finite size effects with respect to isospin ones may be 
estimated considering the caloric curves with the same value of A corresponding to the symmetric 
nuclei (100,50) and (300,150). 

of the curves suggests that for most of the presently analyzed multifragmentation reactions 
finite size effects are small enough to be safely negligible. 

The monotonic decrease of the critical temperature and pressure with the isospin asym- 
metry, together with the reduction of the coexistence region, are illustrated in Fig. [21 where 
the phase diagrams of the 200-nucleon systems with different asymmetries are projected in 
the temperature-total energy, pressure-temperature and pressure-total energy planes. The 
solid lines correspond to the borders of the coexistence region, obtained from a Maxwell 
construction on the constant A microcanonical caloric curves. The dashed lines indicate 
the borders of the spinodal zone, defined by the back-bending extension for each A value. 
The thick solid line connects the critical points in any representation. By extrapolating 
this line, we can see that pure neutron and proton systems may exist in the super-critical 
phase only, as expected from their inability to form clusters at any temperature. It is inter- 
esting to remark that this intuitive result is obtained only as a limiting situation, while a 
phase transition survives with a sizeable critical temperature for systems as asymmetric as 
Z/A = 0.25. 

These results are in qualitative agreement with nuclear matter calculations 0, I^, Q] , show- 
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FIG. 2: (Color online) Projections in the temperature-total energy, pressure-temperature and 
pressure-total energy planes of the phase diagrams of 200-nucleon systems with various isospin 
asymmetries {6 = (N — Z)/A=0, 0.18, 0.30 and 0.50). The solid lines mark the borders of the 
coexistence region while the dashed lines indicate the borders of the spinodal zone. The dotted 
lines in the upper panel show the paths followed through the phase diagram by the considered 
systems when the average freeze-out volume is fixed to 6Vo and the excitation energy ranges from 
2 to 10 MeV/nucleon. The full symbols point the location of multifragmentation events with £'ex=6 
MeV/nucleon further discussed in the text. 
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ing that the thermodynamics of fragmentation of a finite nuclear system can be associated 
to the phenomenology of the nuclear matter liquid-gas phase transition. This intuitive con- 
nection is systematically pushed forward in experimental studies, however from a theoretical 
point of view this is not a trivial issue. If caloric curves and heat capacities in the Statis- 
tical Multifragmentation Model (SMM) are available since more than two decades, most 
calculations with finite systems 38| were performed at constant volume and not constant 
pressure, and therefore lead to signals that cannot univocally be interpreted as a first order 
phase transition. A detailed exploration of the fragmentation phase diagram was presented 
in Refs. jsl, . The results show a continuous transition from a high temperature single 
phase to a mixed phase, this latter extending over the whole density domain of validity of 
the model This is reminiscent of the liquid-gas phenomenology, but the transition from 
the liquid side, which is the only one accessible experimentally, cannot be studied within 
this analytical model, and is hereby presented for the first time. The isospin dependence of 
the fragmentation phase diagram has never been studied before to our knowledge either. 

IV. ISOTOPIC DISTRIBUTIONS AND FRACTIONATION 

In this section we turn to the connection between the system phase diagram and cluster 
observables. 

The generic feature of a first order fluid phase transition with two conserved particle 
numbers is the fractionation phenomenon: if the coupling between alike particles is less 
attractive then the coupling between not alike ones, the ordered phase is systematically 
more symmetric than the disordered one. In the previous section we have shown that the 
fragmentation transition is qualitatively similar to such a fluid transition, we can therefore 
expect to find a trace of fractionation in the fragments and particles chemical composition. 

Looking at the isotopic composition of fragments of different size emitted by neutron 
rich nuclei, it has been observed that the isospin ratio A/Z is a monotonically decreasing 
function of the fragment size, as one would intuitively expect if fractionation takes place 
and thermodynamic discontinuities are rounded by finite size effects. This "fractionation" 



phenomenon has been observed in the experimental analyses 18|, ll9|, l20|, l2l|, |22] , but also 



in different dynamical models 15|, ll6|, ll7| , where it does not seem always connected to the 



phase coexistence phenomenology. 
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A = 200, <V> = 6 Vo, E„ 
A = 6 A----. 
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FIG. 3: (Color online) Break-up isotopic yield distributions of isobars with A=6, 15, 20 and 30 
originating from the multifragmentation of 200-nucleon sources with different isospin asymmetries 
at 6 MeV/nucleon excitation energy and an average freeze-out volume {V) = 6Vq. 

As we have stressed in section HI when dealing with finite systems the properties of 
coexisting phases cannot be deduced from the properties of pure phases by a simple linear 
combination . It is therefore not clear whether the neutron (proton) enrichment of the nuclear 
gas (liquid) characteristic of phase coexistence in neutron rich nuclear matter ^, Q] will 
be apparent in the partitions of the finite system inside the coexistence region. 

Fig. [3] presents isotopic yield distributions of isobars with A = 6,15,20,30 obtained in 
the multifragmentation of nuclear systems (200,100) {6 = {N-Z)/A=0), (200,82) (5=0.18), 
(200,70) (5=0.30), (200,50) (5=0.50), (200,130) (5=0.30) and (200,150) (5=0.50) in a state 
representative for most multifragmentation reactions, Eex=Q MeV/nucleon and {V) = GVq, 
where Vq is the volume corresponding to normal nuclear density, while the free neutron and 
proton multiplicities are listed in Table IIV[ These states are located in very different regions 
of the phase diagram. This is shown in the upper panel of Fig. [2l where the dotted lines 
mark the paths followed by the considered systems when the average volume is fixed to 
(V) = 6Vo and the excitation energy increases from 2 to 10 MeV/nucleon. We can see that 
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Multiplicity / Source 


(200,100) (200,82) 


(200,70) 


(200,50) 


(200,130) 


(200,150) 


neutron 


6.71 


1.31 -10^ 


2.09 -10^ 


3.27 -10^ 


1.28 


1.80 -10-1 


proton 


6.71 


3.10 


1.27 


1.80 -10-1 


2.07 -10^ 


3.24 -10^ 



TABLE I: Neutron and proton break-up multiplicity for different 200-nucleon systems with an 
excitation energy of 6 MeV/nucleon and an average freeze-out volume {V) = 6Vq. 

the {V) = GVq and Eex=6 MeV/nucleon state (reported by full circles) is situated well inside 
the spinodal zone for 6=0 and 5=0.18, while it is close to the critical point for 5=0.30 and 
belongs to the super-critical region for 5=0.50. 

The distributions of Fig. [3] exhibit some trivial characteristic features: the isospin sym- 
metric source produces preferentially isospin symmetric break-up fragments (the isotopic 
yield distributions have a maximum at Z = A/2), and break-up fragment formation is in- 
variant to n-p inversion (neutron and proton yields are equal and isotopic yield distributions 
are symmetric with respect to A/2). Concerning the isospin asymmetric systems, one can 
see that the more neutron (proton) rich is the source, more free neutrons (protons) are 
emitted and neutron (proton) richer are the break-up fragments. The isospin invariance in 
the absence of Coulomb is confirmed by the isotopic yield distribution of the mirror nuclei 
(200,70) vs. (200,130) and (200,50) vs. (200,150), which have reflection symmetry with 
respect to the Z = A/2 axis. 

One can also notice that the distributions are cut both on the proton rich and on the 
neutron rich side. This is due to the dramatic decrease of the binding energy with isospin 
asymmetry approaching the drip-lines. As a consequence, for asymmetric sources primary 
fragments tend to be more symmetric than the initial source, the total asymmetry of the 
system being preserved by a correspondingly increased number of free neutrons (protons). 
The fractionation induced by this effect can be appreciated from Fig. HI which shows as a 
function of the fragment size its average isospin content for the four different asymmetries 
considered above. 

The first feature arising from Fig. H] is that fragments are usually more proton-rich than 
the corresponding source. The only exception corresponds to the symmetric source and 
fragments whose charge is close to half the source charge, where mass and charge conservation 
induces for the Z/A ratios values slightly lower than 0.5. Thus, the approximation frequently 
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invoked in multifragmentation studies, that primary fragments have the same NjZ of their 



emitting source 



22 



25|, 



46 



47l |. does not seem to be correct, if primary partitions correspond 



to statistical equihbrium. 

Even more important, the degree of fractionation is seen to monotonically increase with 
the asymmetry of the source, independent of the location of the multifragmentation event in 
the phase diagram. Indeed, the occurence of fractionation directly follows from the isospin 
content of the free particles because of mass and charge conservation and, thus, it cannot 
be taken as a signature of coexistence in finite systems. On the other hand, the behavior 
of fractionation with the fragment size is very different depending on the thermodynamic 
characterization of the system. Inside the spinodal region (two upper panels), the (ZIA) vs. 



Z distributions shows a clear " U" shape which has been already discussed in Ref . [48|] . The 
two bottom panels which refer to system in a "pure" phase (liquid or fluid) present similar 
characteristics which are very different from the behavior discussed above: both show a 
monotonically increasing {Z/A) vs. Z distribution which for Z source/^ reaches saturation at 
about {Z/A) ^0.42. 

These observations mean that the fractionation phenomenon naturally appears as soon 
as the fragmentation process is ruled by thermal laws. It allows to identify the coexistence 
region of the first order phase transition only if an accurate isotopic characterization of all 
emitted fragments is possible. 

The energy and asymmetry dependence of fractionation is further explored in Fig. [5l 
which gives the evolution with excitation energy of various ratios of light mirror nuclei 
isotopic yields. 

The top panel corresponds to the symmetric source (200,100) and the results indicate 
that, no matter the excitation energy, mirror nuclei are produced with equal probability. 
The lower panels show that this is not true in asymmetric systems, and the increase in 
the emission probability for asymmetric light clusters respect to combinatorial expectations, 
increases with the asymmetry of the source. Similar to the results of Fig. HI we can see 
that fractionation is mainly dictated by the number of evaporated nucleon in excess, and 
no special pattern may be distinguished for the events located inside the coexistence region 
respect to those situated in the liquid or supercritical regimes. Indeed, similar values are 
obtained, for instance, for the 5 = 0.18 source with Ef.x=4: MeV/nucleon (phase coexistence 
region) and the 6 = 0.30 source at Eex=8 MeV/nucleon (supercritical region). Concerning 
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FIG. 4: (Color online) Average isospin content (Z/A) of the break-up fragments produced in 
the multifragmentation of different neutron rich 200-nucleon systems with asymmetries 6=0, 0.18, 
0.30, 0.50 as a function of fragment charge. In all cases (V) = 6Vo and Eex=Q MeV/nucleon. The 
horizontal solid lines indicate the isospin content of the sources. 

the dependence with excitation energy, we can see that the energy increase partially washes 
out the trend of the neutron rich systems to preferentiall y p roduce neutron rich fragments. 



This is in agreement with statistical microscopic models [49 1. 

Based on the idea of isospin fractionation, it has been proposed that the gas neutron 
enrichment can be measured from such ratios 
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25|. Indeed, in the grand-canonical ap- 



proximation, if the charge difference between the two isobars is AZ = 1, the isobaric ratio 
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FIG. 5: (Color online) Ratios of isotopic yields of different light mirror nuclei at the break-up stage 
of 200-nucleon systems with different isospin compositions as a function of excitation energy. In 
all cases the average freeze-out volume is {V) = GVq. For the neutron-rich sources the considered 
ratios are with the neutron-rich isobars in the numerator while for the neutron-deficient sources 
the neutron-rich isobars are in the denominator. 
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is given by 

where AB = B{A, Zi) — B{A, Z2), and A/x = /i„ — /ip, if the neutron rich isobar is in the 
numerator. Choosing hght isobars having close binding energies, this ratio is then a direct 
measure of the chemical potential difference between neutrons and protons, i.e. of the ratio 
Pn/ Pp of the free neutron-proton densities. 

Fig. shows that, at least for the considered fragments, this approximation is reasonable 
enough for all excitation energies and asymmetries, and information on the free neutron 
versus proton behavior can indeed be inferred from the measurements of isotopically resolved 
light fragments, assuming that freeze-out yields can be restored from the experimentally 
detected cold fragments. However, it is important to stress that this free densities ratio 
cannot unambiguously sign the neutron enrichment of the gas phase, since, as discussed 
above, the same behavior is observed when no gas phase can be thermodynamically defined. 

To conclude, in this section we have pointed out that isospin fractionation cannot be taken 
as a signature of phase coexistence when dealing with finite systems. In a previous work 



2J], we have already shown that the Coulomb interaction tends to quench the coexistence 
zone: because of that the multi-fragmentation phenomenology can be associated with a 
super-critical region of the charged-system phase diagram. Here we show that, even in the 
absence of the Coulomb interaction, a generic universal feature of fragmentation, namely 
isospin fractionation, can show- up above the critical point. If the effect of Coulomb on 



the phase diagram strongly depends on the specific model used [2J, |52| , this super-critical 
fractionation is a generic effect, which we believe should be present in any fragmentation 
model. Indeed it is due to the combined effect of clustering in the super-critical region (which 
favors the formation of fragments close to stability, i.e. an isospin symmetric fraction of the 
system condensed at finite baryon density), and particle number conservation (which forces 
the low density non-clustered part to have a strong isospin asymmetry). These features are 
naturally present in any model, microscopic or macroscopic, respecting conservation laws 
and ruled by the competition between entropy and energy. Such an effect is not accessible 
in nuclear matter calculations, where the sharing of the system into a dense and a diluted 



fraction is by construction a sign of phase separation. Classical models 53|] have already 
shown clustering in the supercritical region, but this is to our knowledge the first time that 
this effect is reported in a realistic nuclear multifragmentation model. 
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V. ISOTOPIC WIDTHS AND SYMMETRY ENERGY IN FINITE NEUTRAL 
SYSTEMS 



A very powerful motivation in the study of isotopic distributions in fragmentation reac- 
tions is given by the well-spread expectation that information coming from the low density 
finite temperature coexistence zone of the phase diagram will be sensitive to the symmetry 
energy coefficient of the nuclear (free) energy - de nsity functional at finite temperatures, 



and densities well below saturation 
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261, 



23, 



28l . l29l . l50|]. Such analyses would then be 



complementary to isospin diffusion and neutron skin measurements 5l|] and additionally 
give a unique information on temperature effects on the symmetry energy. 

In the MMM model the binding energy of a cluster of mass A and charge Z is parame- 
terized as, 

{A - 2Zf 



Bnoc{A,Z) = (a^A-asA^/^) -a,{a^A-asA^/^)- 



A^ 



(a„A 



^A ^'^) — Csym{A) 



(A - 2zy 

A 



(3) 



and includes a full mass dependence of the bulk+surface and isospin-dependent contribu- 
tions. The Coulomb part of the binding energy {acZ'^ / A"^ + aaZ"^ / A) 5^] which is included in 
the standard version of MMM 23(], is switched off for this study, to concentrate on isospin 
effects. 

As we have stressed in section [Til in the framework of statistical models under the Fisher 
approximation 35|], the low density correlations are entirely exhausted by clusterization. 



This means that the symmetry (free) energy entering in the fragment production yields 
inside coexistence should be the symmetry energy of isolated nuclei at finite temperatures. 
In particular the interaction part of this energy (Eq. ([3])) should correspond to normal ground 
state values, meaning that the liquid drop parameters ai,ay,as have standard ground state 
values [5J]. In our model isotopic yields are therefore expected to be entirely determined in 



the whole phase diagram by the symmetry energy coefficients Csym{A). 

If this approximation gives a correct description of multifragmentation, this would mean 
that no relevant information on Csym{p, T) can be inferred from fragment observables. If, on 
the other hand, the energy functional of break-up fragments differs from the ground-state 
functional 26|, |27|, it would be interesting to trace its behavior and isotopic distributions 



could be good candidates. The extent to which this may be true is a difficult theoretical 
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issue, demanding a quantal many-body transport treatment, completely out of scope of 
the present study. Whatever the final answer, for the fragment symmetry-energy to be 
accessible from experimental data, it is necessary to prove that in a controlled model where 
the symmetry energy set in break-up fragments is an input value of the calculation, the 
proposed observables do indeed recover its value within a good precision. 

The information on Csym can be directly inferred from the widths of the isotopic dis- 
tributions in the grandcanonical approximation. Indeed, a Gaussian approximation on the 
grandcanonical expression, 

>^A,„,,,(iV, Z) = Z-^^^^ exp [-(3 (F^(iV, Z) - /.„iV - , (4) 

(N-Z-Io) 



leads to 



Yp,,.{A,N-Z) = K{A)exp 



(5) 



2aj{A) 

where Jq = — Z is the most probable value of N — Z for a given value of the cluster size 
A, K{A) does not depend on the asymmetry I = N — Z and the isospin variance is related 
to the symmetry energy coefficient by, 

ajiA) = . (6) 

The coefficient Cfj^^(A) appearing in this last expression is a free symmetry energy coef- 
ficient given by, 

Cfym(^) = U, (7) 

and coincides with Csym{A) defined by Eq. ([3]) if we neglect the I dependence of the 
excitation energy and entropy associated to a given mass A. In this case Eq. [6] reads, 

2 

"^CsymiA) 

The quality of all these approximations can be appreciated from Fig. [6l which displays the 
width of the asymmetry distribution (open circles) as a function of the fragment mass for 
the same thermodynamic conditions as Figs. El IH El For all considered asymmetries, Eq. 

(solid line) appears well verified for small masses, meaning that this width can indeed be 
taken as a measure of the underlying symmetry energy. For higher masses and/or extreme 
asymmetries the Gaussian approximation breaks down as it can be observed from Fig. [3l 
and the link between sjTiimetry energy and fiuctuation is lost. The neutron distribution 
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(filled circles) contains approximately the same information brought by the asymmetry dis- 
tribution. The dashed lines in Fig. [6] give the grandcanonical expectation for the isotopic 
widths Eq. (IHD, when only the bulk term (C^^^ = aja^,) of the symmetry energy is consid- 
ered, instead of the complete expression {CgymiA) = ata^ — ajOsA^^/^) used for the solid 
line. We can see that accounting for surface effects in the fragment energy functional leads 
to a considerable increase of the isospin widths even for relatively massive fragments. Such 
fragments {A > 30) are not adapted to study the symmetry energy coefficient though, be- 
cause of the important effect of the mass and charge conservation constraint that causes the 
widths to deviate from their grandcanonical expectation. 

It is very interesting to observe that the solid and dashed lines are almost parallel to 
each other. This means that the expected functional dependence of the width on the mass 
number does not change drastically if the symmetry energy has a surface contribution or 
not. As a consequence, a large width as expected for light, surface-dominated fragments 
{A < 30), can be easily mis- interpreted as a signature of a reduced bulk symmetry value, as 
it has suggested by some recent publications Q ■ 

The relation between isotopic widths and symmetry energy has interesting consequences 



on the isoscaling observable, which has raised a great interest in the recent literature 



a 
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5o| . As long as the distributions can be approximated by Gaussians, the ratio between 
the production yield of the same isotope in two different systems (1) and (2) (where we 
denote by (2) the neutron-rich one) for a given Z can be expressed as a function of neutron 
number N as: 



Yii)iN,Z)J 2 



where A^(j) is the most probable N value for the element Z in system (i) and a^^.^ is the 
variance of the distribution in the same system. If we choose as systems (1) and (2) 
two systems of similar masses and temperatures, then aj^ ~ aj^ ~ cr^ and the ratio 
shows (in log scale) a linear dependence on A^ (at fixed Z). Similar arguments hold also for 
fragments with fixed A^, meaning that the quantity in left-hand side of Eq. (Q has also a 
linear dependence on Z (at fixed A^). This result is known in the literature as the isoscaling 



phenomenon 251 ]: 
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FIG. 6: (Color online) Widths of Y{I')\a (open symbols) and Y{N~)\z (full symbols) distributions 
corresponding to the break-up stage of 200-nucleon systems with different asymmetries (5=0, 0.18, 
0.30, 0.50) as a function of fragment mass in comparison with predictions of Eq. ([8|) calculated for 
fragment (full line) and bulk (dashed line) symmetry energies. All systems are characterized by 
(y) = 6Vo and i?ex=6 MeV/nucleon and the Coulomb interaction is switched off. 
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where, 

a(Z) = -l(iV(2)-iV(i)), (11) 

or, for symmetric distributions, 

a{Z) = -1 ((iV)(2) - (iV)(i)) . (12) 

In the mass region where the Gaussian approximation is well verified and in the absence 
of Coulomb effects (see Fig. [H]), 

aUZ)^a]{{A){Z)). (13) 
Then, the isoscaling parameter a{Z) is linked to the fragment symmetry energy by, 

c^iZ) ^ ^^^Ml ((iv)(,) - (iV)(,)) , (14) 



or, equivalent ly. 



ai< Z{A) >) ^ '^.£1^ ((/)(,) _ (/)(,)) . (15) 



Fig. [7] compares the symmetry energy coefficient extracted from Eq. (IMI) with the 
input symmetry energy of the model for a representative case. Different average volumes, 
excitation energies, asymmetry ratios give similar results. We can see that once again the 
Gaussian approximation appears well verified for light fragments and, in that case, isoscaling 
technics give a reasonably good measure of the fragment symmetry energy, the reason being 
that this parameter is directly linked to isotopic fiuctuations. 

A similar expression, 

, , ^ ic.,„m) ( 2^ (iji 



has been derived in Ref. 16j from Eq. (jl]) with a similar saddle point approximation as 
for Eq. (IT^ . considering only the most probable isotopes for each Z. This equation is also 
plotted on Fig. [7] and gives comparable results to Eq. f|T^ . 

In the some experimental analyses the isoscaling parameter a{Z) is not extracted sepa- 
rately for each isotope, as defined in Eq. ( fTOj) . but as the average value over fragments with 
1 < Z < 8. For this last situation we adopt the notation aexp- Fig. [7] plots the behavior of 
Eq. (fT6|) for these different definitions of a. It comes out that the procedure to calculate 
a using exclusively the light nuclei does not perturb the extracted symmetry energy, but 
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FIG. 7: (Color online) The input symmetry energy coefficient as a function of fragment size (solid 
line) is compared to the estimations from fragment observables, calculated for the break-up stage of 
(200,82) and (200,100) nuclear systems without Coulomb with {V) = 6Vo and £'ea;=6 MeV/nucleon. 
Open diamonds and circles stand for predictions of Eq. ([8]) and, respectively, Eq. (jl4p . Predictions 
of Eq. ()16p (stars) are plotted as open or full stars depending whether the isoscaling parameter 
a is calculated according to its definition (Eq. [TO]) or as the average values of fragments with 
1 < Z < 8. 

rather minimizes the deformations due to conservation laws, which would prevent a precise 
extraction of the symmetry energy coefficient for large clusters. Quantitatively speaking, by 
ignoring the monotonic increase with Z, the use of aexp results in slightly lower values of 
the symmetry energy with respect to the ones obtained employing a{Z). 

It has been recently argued 0, [s^ that a constant value of the isoscaling parameter 
a would imply a bulk character for the associated symmetry energy. It is particularly 
interesting to notice that in our model this is not the case. Indeed, the size dependence of 
Csym{{A))/ (A) is compensated by the size dependence of ((A^)(2) — (^)(i)) giving a constant 
a. 
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VI. ISOTOPIC WIDTHS AND SYMMETRY ENERGY IN REAL NUCLEI 



The robustness of these signals to measure the symmetry energy when the Coulomb 
interaction is included is particularly important as it gives the extent to which one may 
extract this basic quantity from multifragmentation data, assuming that one may access the 
chemical composition of the physically relevant break-up fragments. 

In ground state nuclei the Coulomb interaction is known to shift the stability peak to- 
wards neutron rich nuclei and to reduce dramatically the binding energy of the proton rich 
ones. This last effect is responsible for a strong narrowing of the B{A, Z)\a and B{A, Z)\isr 
distributions and gets more pronounced with the mass increase. Since the logarithm of 
fragments multiplicity approximately follows the evolution of B{A, Z)/T (Eq. (jlj)), we can 
expect to find the same effect of a width reduction on fragments yields. 

Fig. [Hlplots the widths of the Y{I)\a, ^(A^)!^ and, for the sake of completeness, Y{Z)\jq 
distributions for two nuclei, (210, 82) and (190, 82), in a thermodynamical state relevant 
for most multifragmentation reactions, V = 4Vo and Eex=Q MeV/nucleon. The choice of 
a pair of nuclei with equal charge minimizes the interference between Coulomb and isospin 
contributions. 

We can observe that the expected dispersion between cr|, a% and cr| increases with 
the fragment mass and the source asymmetry, which obviously restricts the validity of the 
approximation Eq. (fT3l) . The decrease of a] and (t| can be attributed to the above mentioned 
narrowing of B{A,Z)\a and B{A,Z)\n distributions under Coulomb effect. The increase 
of aff then results from particle number conservation. Indeed, under mass and charge 
conservation the dispersion of fragment yields obeys the law: 

^ {al + 4) /2. (17) 

The consequence of these effects is that the slight deviation of aj from Eq. ([H]), already 
present for the heavy fragments obtained in the decay of finite neutral systems, gets more 
pronounced under Coulomb such that for A > 50, a] practically falls over the predictions 
of Eq. ([8]) with bulk symmetry energy (dashed line) instead of fragments symmetry energy 
(solid line). However it is important to observe that, for the light fragments which are 
isotopically resolved in most experimental data-sets, aj stays an excellent measure of the 
input symmetry energy irrespective of conservation laws and Coulomb effects. 
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FIG. 8: (Color online) Widths of Y{I)\a (open circles), Y{N)\z (full circles) and Y{Z)\n (open 
stars) distributions corresponding to the break-up stage of two Z=82 nuclei and different asymme- 
tries, (5=0.14 (upper panel) and 0.22 (lower panel) characterized by the freeze-out volume V = 4Vq 
and excitation energy E'ex—S MeV/nucleon as a function of fragment mass in comparison with 
predictions of Eq. ([8|) calculated for fragment (full line) and bulk (dashed line) symmetry energies. 

Finally, Fig. Ogives the quality of the approximations of the different formulas (Eqs. ([HD, 
ffn|) . ffT^ and ffTBl) ) proposed in the previous section to access the symmetry energy. As 
in the previous section, for the so-far widely used Eq. f|T6l) we considered two 'definitions' 
of a: the slope of the logarithm of isotopic yield ratios for each Z (Eq. (ITO!) ) and, in 
the spirit of most experimental analyses, the average value of a{Z) for fragments with 
1 < Z < 8. Not surprisingly, the reconstruction of the input symmetry energy is perturbed 
by Coulomb effects. As one may notice, for heavy fragments, the violation of the different 
approximations on which the proposed expressions reckon leads to a relative dispersion 
which may range from 15% (for A = 20) to 25% (for A = 70). The best description of the 
fragment symmetry energy seems to be offered by the simplified version of Eq. f|T5l) where 
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mass 

FIG. 9: (Color online) The same as in Fig. [3 but for the break-up stage of (210,82) and (190,82) 
nuclei with V = 4Vo and Eex=(^ MeV/nucleon when the Coulomb energy is included. The esti- 
mation of Eq. (open squares) is also presented. For this last expression, the abscissa has the 
meaning of the exact mass. 

the parameter a is calculated as the average value of Eq. ffTOl) over fragments with Z=l 
to 8 as it underestimates the witness value by less than 1 MeV over the whole considered 
fragment mass interval [A = 1 — 70). As a general statement, considering the light {A < 20) 
fragments used in most experimental analyses for which the conservation laws constraint are 
the least important, it is encouraging to see that both isoscaling and isotopic widths give 
an consistent estimation of Csym which deviates from the input value of no more than 20% 
55|. 

VII. CONCLUSIONS 

Isospin effects on the thermal and phase properties of finite uncharged excited nuclear 
systems at sub-nuclear densities have been studied in the framework of a microcanonical 
statistical model with cluster degrees of freedom. We find that the isospin asymmetry of the 
source reduces the width of the coexistence region and the critical temperature and pressure 
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values, in qualitative agreement with the well known results for infinite nuclear matter. The 
similarity of the phase diagram with the nuclear matter one is an extra confirmation that 
the multifragment production phenomenon can be associated to the coexistence zone of a 
first order phase transition of the liquid-gas type. 

To push this connection further, we have explored the relation between fragment chemical 
composition and the expected isospin fractionation in the coexistence region of a multi-fiuid 
system. We have shown that a number of excess neutrons are emitted as free particles in 
neutron-rich sources. This number strongly increases with the asymmetry of the source, 
independent of the system location in the phase diagram. This implies that free nucleons 
emitted by a finite isolated system cannot be unambiguously associated to the gas phase 
of the corresponding phase diagram. Because of the mass and charge conservation law, 
the isotopic distribution of complex fragments is in turn mainly dictated by the excess 
free nucleons. Thus the fractionation phenomenon cannot be taken as a measure of phase 
coexistence. An interesting observable is given by the 'U-shape' of the complete ZjA vs. Z 
distribution which appears characteristic of the phase coexistence region, while more simple 
quantities like ratios of mirror nuclei give ambiguous results. 

Finally, wc investigate the relation between fragments symmetry energy and the variance 
of isotopic distributions. A simple expression relates the symmetry energy with the isospin 
asymmetry variance and a new formula is proposed in order to calculate the same quantity 
out of isoscahng observables. While they are always accurate enough in the case of finite 
neutral systems, some of these expressions show sizable deviations in real nuclei under 
the infiuence of Coulomb. The present study, however, suggests that the widths of the 
isotopic distributions, as well as the isoscaling parameter, can still give a correct estimation 
of the symmetry energy in physical nuclear multifragmentation data, provided the break-up 
fragment partitions can be restored from the detected cold fragments. 
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